Angular rigidity in tetrahedral network glasses 
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A set of oxide and chalcogenide tetrahedral glasses are investigated using molecular dynamics 
simulations. It is shown that unlike stoichiometric selenides such as GeSe2 and SiSe2, germania 
and silica display large standard deviations in the associated bond angle distributions. Within 
bond-bending constraints theory, this pattern can be interpreted as a manifestation of broken (i.e. 
ineffective) oxygen bond-bending constraints. The same analysis reveals that the changes in the Ge 
composition affects mostly bending around germanium in binary Ge-Se systems, leaving Se-centred 
bending almost unchanged. In contrast, the corresponding Se twisting (quantified by the dihedral 
angle) depends on the Ge composition and is reduced when the system becomes rigid. Our results 
establishes the atomic-scale foundations of the phenomelogical rigidity theory, thereby profoundly 
extending its significance and impact on the structural description of network glasses. 
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The large variety of physico-chemical behaviours in- 
herent in tetrahedral network glasses (in particular, those 
involving Group IV (A=Si, Ge) oxides (X=0) or chalco- 
genides (X=S,Se,Te)) is deeply related to the underly- 
ing network topology, i.e. the nature of the connec- 
tions (edge- or corner-sharing) among the basic tetra- 
hedral structural units [3|, |4|. In the search of an uni- 
fying approach, it is tempting to follow rigidity theory, 
that describes the interplay between network properties 
and connectivity by considering covalent networks in very 
much the same fashion as mechanical trusses [5|- This is 
achieved via enumeration of mechanical rigid constraints 
nc arising from relevant atomic interactions, r/2 bond- 
stretching (radial) and (2r — 3) bond-bending (angular) 
for a r-folded atom. As a consequence, a rigidity transi- 
tion at the network mean coordination number f=2.4 is 
identified '5y , separating underconstrained networks hav- 
ing low-frequency (floppy) deformation modes (uc < 3) 
from overconstrained ones (uc > 3). 

Accordingly, all stoichiometric compounds AX2 should 
be stressed rigid i.e. they have more constraints than 
degrees of freedom (ric— 3.67), with a mean coordina- 
tion number of f=2.67 larger than the critical fc—2A 
for which optimal glass formation with small enthalpic 
changes at the glass transition, is supposed to occur [6|. 

The legitimacy of such a picture is challenged by the 
observation that silica (Si02) and germania (Ge02) are 
found to form rather easily glasses [7|, in contrast with 
the corresponding chalcogenides (e.g. GeSe2) found at 
the very limit of the binary glass-forming region (in e.g. 
Ge^Sei-a; [8|). Furthermore, oxides have a low frequency 
(floppy) contribution in the vibrational density of states, 
suggesting that these systems are flexible [9J or, at least, 
nearly optimally constrained [10| (isostatic, nc=3). The 



latter result can be recovered in rigidity theory under the 
heuristic assumption that enhanced oxygen bond-angle 
values lead to broken angular constraints |ll| reducing Uc 
from 3.67 to 3.0. However, a microscopic scale rationale 
for this specific network behavior is lacking. 

These pieces of evidence call for the basic assertions 
of rigidity theory, when applied to the above oxide and 
chalcogenide networks, to be revisited. Here we rely on 
an atomic-scale approach (as molecular dynamics) able 
to substantiate and enrich the general trends of rigidity 
theory via the explicit account of the details of chemical 
bonding. 

In this Letter we bridge the gap between constraint 
counting algorithms and the statistical mechanical be- 
havior of relevant atomic-scale quantities (angular dis- 
tributions) by considering a set of chalcogenides and ox- 
ide glasses. Two issues commonly tackled within rigidity 
theory are given new atomic-scale foundations. First, 
by focussing on the stoichiometric AX2 composition (i.e. 
Si02, Gc02, SiSe2 and GcSe2) we show that oxide sys- 
tems contain tetrahedra which act as rigid units having 
a much smaller minimal angular excursion (quantified 
by ae , the standard deviation of the partial bond angle 
distributions) than their chalcogenide counterparts. By 
analysing in a similar fashion the binary Gea;Sei_a; sys- 
tem with changing compositions, we find that changes 
in rigidity (increasing Ge content x) are noticeable in 
the angular environment of germanium, leaving the Se 
centred angular excursions nearly unchanged when mov- 
ing from a flexible to a stressed rigid phase. It appears 
that the amplitude of the angular distribution around A 
atoms is increased in the stressed rigid phase with an 
increased distorsion (i.e. higher cr^s) of the tetrahedra. 
This result allows reconsidering the general accepted pic- 
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FIG. 1: (color online) From top to bottom: oxygen, selenium 
and germanium partial bond angle distributions (PBAD) in 
Ge02 and GeSe2 for an arbitrary N=6. The colored curves 
correspond to PBADs having the lowest standard deviation(s) 
ae (see Fig. 2). The sharp peaks at ~ 40° correspond to 
the hard-core repulsion. Labels defined in the bottom panel 
are used throughout the text. 



ture of a Ge-Se network made of flexible Se-chains and 
rigid GeSe4/2 tetrahedra [5|. As a second main outcome, 
we find that ag for oxygen in oxides, is much larger than 
ag for selenium in GeSe2, proving that oxygen bending 
constraints are broken (i.e. ineffective) in Si02 and Ge02. 
This provides a microscopic rationale for the sensitivity 
to bond-bending around oxygen in oxide networks, con- 
sistent with experimental evidence 
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Our analysis is based on molecular dynamics trajecto- 
ries obtained at T=300 K for a set of glassy systems 
encompassing Ge02, Si02, SiSe2 and Ge^Sei-^; for 4 
different compositions: x = 0.10 (GeSeg in the flexi- 
ble phase), X = 0.20 (GeSe4 at the rigidity transition 
[6|), X = 0.25 (GeSca in the intermediate phase [8|) and 
X = 0.33 (GeSe2 in the stressed rigid phase). Oxides at 
the experimental densities have been simulated accord- 
ing to Ref. [12| and 13 1 using a classical Born-Mayer 



force field. Due to the large difference of electronega- 
tivity between Si(Ge) and O, this choice ensures plausi- 
ble qualitative modelling within classical molecular dy- 
namics. In the case of SiSe2 and Ge^Sei-^; systems, 
we resort to first-principle molecular dynamics (FPMD) 



FIG. 2: (color online) Standard deviation ag. . of partial 
bond angle distributions (PBADs) in oxides (Si02, Ge02, top 
panel) and chalcogenides (SiSe2, GeSe2, bottom panel) as a 
function of the angle number (see defintion in Fig. 1). 



within a fully self-consistent framework that proved ade- 
quate to describe chemical bonding and its changes with 
concentration. For some of the systems, temporal tra- 
jectories recorded previously are substantially extended 
to attain optimal statistical accuracy (within a few per- 
cent at most). Otherwise, as in the case of GeSeg and 
GeSca, glassy structures are produced from the outset 
after cooling from the liquid state and appropriate struc- 
tural relaxation. Overall, typical time trajectories for col- 
lection of the averages in the glassy state cover ^ 100 ps. 
Details on the FPMD methodology and the productions 
of g lassy structure are extensively reported in [iJ] and 
[l5(. For all glassy structures, the atomic-scale picture 
is consistent with that obtained by using an alternative 
first-principles approach [l6[. 

To obtain the number of bond-stretching interactions 
we have focused on neighbour distribution functions 
(NDFs). A set of NDFs can be defined by fixing the 
neighbor number n (first, second etc), the sum of all 



NDFs yielding the standard pair distribution function 
gi(r) (i=A,X). Integration of gi(r)(i=A,X) up to the first 
minimum gives the coordination numbers vx and r^, and 
hence the corresponding number of bond-stretching con- 
straints ri/2 {i—A,X). For all considered systems, we find 
rx — 2 and rA = 4: leading to 1 and 2 bond-stretching 
constraints for the X and A atoms. For instance, when 
integrating the pair distribution function up to its first 
minimum, we find for Ge02 rce—^-Ol- and ro = 1.97, and 
for GeSe2, rQf=4:.02 and rse=l-96, in agreement with 
experiments [4l.ll7l.ll8|. 

Bond-bending constraint counting is based on partial 
bond angle distributions (PBADs) P{Oij) defined as fol- 
lows: for each type of central atom 0, the N first neigh- 
bours i are selected and the N(N-l)/2 corresponding an- 
gles iOj (i=l..N-l, j=2..N) computed, i.e. 102, 103, 
203, etc. The standard deviation ag^- of the distribu- 
tion T'{9ij) gives a quantitative estimate about the angu- 
lar excursion around a mean angular value, and provides 
a measure of the bond-bending strength (small values, 
i.e. intact bond-bending constraint) or the bond-bending 
weakness (large values, ineffective or broken constraint). 

Figure 1 shows such PBADs for oxygen (Ge-O-Ge an- 
gle) and germanium (0-Ge-O) in Ge02 and for selenium 
(Ge-Se-Ge) in GeSe2 for an arbitrary N=6 leading to 
15 possible PBAD's. Oxygen displays for the principal 
contribution 102 (angle number 1, in red) a broad distri- 
bution centred around the angle 9 — 135°, corresponding 
to the one defined by the two closest (Ge) neighbours of 
oxygen, in agreement with experiments |19j . Secondary 
distributions 103 and 203 show peaks centred at around 

~ 90° and 75°. In contrast, GeSe2 exhibits a much 
sharper distribution for the same 102 contribution, im- 
plying reduced angular excursions as compared to Ge02. 
The contributions to the bimodal GeSe2 distribution at 
80° and 100° can be respectively assigned to edge-sharing 
and corner-sharing tetrahedra [20|, a feature that is ab- 
sent in oxides. Six germanium centred angles (bottom 
panel) are found to have almost the same distribution in 
Ge02, centred at an angle of 109° typical of the tetrahe- 
dral environment. 

The behavior of the standard deviations for the PBADs 
is shown in Fig. 2, where, for sake of simplicity, each 
PBAD is given a distinct number (see bottom part of Fig. 

1 for the definition). For all chalcogenide or oxide sys- 
tems, the PBADs relative to the Group IV (Si, Ge) atom 
have a low standard deviation ag^ - , of the order of 10-20°, 
for instance crg,^ ~ 7° for the 102'pBAD of Ge02. These 
values are much smaller than those of the other distribu- 
tions (105, 106, etc.), found close at ~40°. In addition to 
very low angular excursions around the tetrahedral angle 
of 109° (Fig. 2 top), oxides feature all (jg^. (i,j<4) nearly 
equal for the six relevant (Ge,Si) distributions. A differ- 
ent situation occurs in stoichiometric chalcogenides (Fig. 

2 bottom, red curve) which exhibit increased bending for 
the angles defining the tetrahedra (angle 3: 104, angle 7: 
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FIG. 3: (color online) Standard deviation ace and ase ex- 
tracted from the partial bond angle distributions (PBAD) for 
four selected compositions in glassy Ge^iSei-i, 




X-A-X 



FIG. 4: Definition of the dihedral angle 5 used in the text. 
Note that the upper black-filled atom can X (rx=2) or A 
(rA=4). 



204, etc.). These results exemplifies the difference in the 
bending nature of the tetrahedra in these two families of 
networks, pointing to a higher rigidity of the tetrahedra 
in oxides, as all anglular excursions are maintained at the 
same low value (typically 7°). 

Angle distributions around the Group VI atoms are 
markedly different. Indeed, the selenium standard de- 
viation (7ei2 (angle number 1) is found to be low (11.7° 
for GeSe2) compared to the corresponding oxide system 
(26.7° for Ge02)- Therefore, the restoring effect associ- 
ated with bending is much lower in Ge02, allowing for 
the identification of a broken angular constraint. One 
should also have in mind that the standard deviations 
found around 20° are close to those found in high tem- 
perature liquids (18-25° at 1373 K for GeSe2), where an- 
gular constraints are assumed to be broken 21| . This re- 



System Si02 Ge02 SiSe2 GeSeg GeSe4 GeSea GeSe2 

5(deg) 30.9 30.4 30.9 35.0 40.3 31.9 30.6 
as{deg) 17.7 17.5 23.9 27.4 25.9 20.6 21.6 



TABLE I: Mean dihedral angle 5 and standard deviation as 
of the dihedral angle distribution for the seven investigated 
systems 



duces the number density of constraints to nc=0.33(|ryi- 
3)+0.67(|rj5f-4)=3.00, i.e. oxides are optimally con- 
strained. Because of its large fraction of ES tctrahcdra 
[15| . the situation of SiSe2 appears to be somewhat in- 
termediate with ag-^^ — 19.6°, resulting from contributions 
arising from both ES and CS tetrahedra which have re- 
spectively ineffective and intact Se bond-bending con- 
straints. 

It is of interest to apply the above rationale to the 
Ge^Sei-x family of systems, since both the elastic nature 
(flexible, rigid) and the connectivity are strongly depen- 
dent on composition. By increasing the Ge concentration 
the angular excursion inside the tetrahedra is seen to in- 
crease (e.g. ace moving up to 20° from less than 10° 
in the n=3 (104) PBAD) while leaving the stiffest an- 
gle (102) constant. On the other hand, bending around 
the Se atoms is nearly unchanged and it does not display 
any noticable change when the system becomes rigid. To 
find a structural parameter pertaining to the Se atoms 
and sensitive to changes in composition, we have to re- 
sort to the dihedral angle S around a Se atom (see def- 
inition in Fig. 4). As shown in Table I, the dihedral 
angular excursion takes a value of as=27A° for the flexi- 
ble GeSeg composition and decreases significantly (down 
to 20°-21°) for the compositions of rigid systems GeSea 
and GeSe2. Therefore, the network adapts to the pre- 
dominant presence of Se atoms both by decreasing the 
angular variability inside the GeSe4/2 tetrahedra, and by 
allowing for enhanced twisting along the Se chains. Fi- 
nally, in view of the behaviour of a a (A=0,Se) from both 
Figs. 2 and 3, it becomes clear that the angular motion 
inside Si02 and Ge02 tetrahedra behaves very similarly 
to flexible GeSeg or optimally constrained GeSe4 and 
GeSea networks. On the other hand, these oxide sys- 
tems definitely contrast with the corresponding stressed 
rigid chalcogenides SiSe2 and GeSe2. 

We showed that structural information gathered from 
molecular dynamics is able to provide an atomic-scale 
counterpart to phenomenological constraint counting 
concepts applied to network glasses. This leads to a 
clear picture of the topological differences between sys- 
tems having the same composition but different chemical 
nature and systems made of the same species but differing 
in composition. The differences between glasses (oxides 



and chalcogenides) of same stoichiometry are rational- 
ized in terms of amplitude of the inter- and intratetra- 
hedral bending angular variations. In oxides, the results 
show increased bending aroud the oxygen atom, consis- 
tent with a direct Maxwell constraint counting and the 
non-stressed rigid nature of these glasses. In systems un- 
dergoing a rigidity transition (e.g. GexSei-x), it is found 
that the angles defining the GeSe4/2 tetrahedron soften 
with decreasing Ge weight while the Se angular bend- 
ing is almost unchanged. Flexibility along the Se chains 
is best accounted for in terms of twisting along the Se 
chains, found to increase with Se content. The present 
approach finally provides a general framework for estab- 
lishng bonding constraints in a neat way via model simu- 
lations, and should be used in the future for establishing 
constraint counting algorithms in more complex glassy 
materials. 
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